Strong coupling expansion in a correlated three-dimensional topological insulator 
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Motivated by recent studies which show that topological phases may emerge in strongly correlated 
electron systems, we theoretically study the strong electron correlation effect in a three-dimensional 
(3D) topological insulator, which effective Hamiltonian can be described by the Wilson fermion. We 
adopt 1/r long-range Coulomb interaction as the interaction between the bulk electrons. Based on 
the U(l) lattice gauge theory, the strong coupling expansion is applied by assuming that the effective 
interaction is strong. It is shown that the topological insulator phase changes to the normal insulator 
phase in the strong coupling limit, namely the value of the chiral condensate is not zero as well as 
in the case of graphene. 

PACS numbers: 71.27.+a, 11. 15. Ha, 11. 15. Me, 03.65.Vf 



I. INTRODUCTION 

Recently discovered topologically nontrivial phases 
have attracted many researchers and offered a new direc- 
tion to modern physics 1,2 . Topologically nontrivial phase 
and trivial phase, in the presence of time-reversal sym- 
metry, are distinguished by the Zi invariant 3,4 . Strong 
spin-orbit coupling is known to be essential to realize 
topological phases, since topological phases originate in 
the parity change in the lowest unoccupied band from 
even to odd induced by spin-orbit coupling. Topolog- 
ical phases are characterized by the gapless edge (sur- 
face) states which are protected by time-reveral symme- 
try. In 3D topological insulators, the surface states are 
described by the two-component massless Dirac fermions. 
The bulk states in such as Bi 2 Se3 are described by the 
four-component anisotropic massive Dirac fermions 5 . It 
is known that the surface states are robust against per- 
turbation and disorder 6,7 . What about against electron 
correlation, i.e. Coulomb interaction? This is a natural 
question, because it has been revealed that strong elec- 
tron correlation is important in many systems and may 
induce novel phenomena. 

A novel Mott-insulating phase was found recently in 
an iridate 8 , a 5c?-electron system, and has gathered much 
attention. Remarkably, the phase is induced by the coop- 
eration of strong spin-orbit coupling and strong electron 
correlation. Evolved by this discovery, many studies have 
been done intensively in systems where both spin-orbit 
coupling and electron correlation exist, for the search for 
novel phases induced by them. Especially, it is of inter- 
est that topological phases such as the quantum spin Hall 
insulator 9 and the Weyl semimetal 10 are predicted in iri- 
dates. These results suggest that topological phases may 
emerge in strongly correlated d-electron systems. Preced- 
ing studies mainly focus on the competition between the 
spin or charge ordered phase and the topological phase 
in Hubbard-like models on honeycomb lattices 11-19 , other 
2D lattices 20-24 and 3D lattices 25-28 . Another study on 
the surface Dirac fermions shows that the Dirac fermions 
become massive with finite correlation strength due to 
the spotaneous magnetization 29 . 



On the other hand, the electron correlation effect in 
graphene, a two-dimensional Dirac fermion system, has 
been studied widely. In graphene in vacuum, the cou- 
pling constant becomes effectively large due to the small 
Fermi velocity. It has been predicted that a finite band 
gap is induced in charge neutral graphene in vacuum. 
In such a case, the strong coupling lattice gauge the- 
ory is applied 30-38 . The chiral condensate is the order 
parameter for the insulator-semimetal transition in the 
lattice gauge theory. It is noteworthy that lattice Monte 
Carlo studies show quantitatively correct critical value of 
the coupling strength below which the system becomes 
gapless 31,32 ' 38 (graphene on a SiC>2 substrate is conduct- 
ing). These results motivated us to do this study. 

In this paper, we focus on the strong electron corre- 
lation effect in a 3D Dirac fermion system on a lattice 
which is a simple model describing a topologically non- 
trivial state. We adopt 1/r long-range Coulomb inter- 
action as an interaction between the bulk electrons, be- 
cause the screening effect in Dirac fermion systems is 
considered to be weak due to the vanishing of the den- 
sity of states. This situation is nothing but what is de- 
scribed by the U(l) lattice gauge theory. Therefore, we 
can perform the strong coupling expansion of the lat- 
tice gauge theory by assuming that the effective coupling 
constant is large. The procedure is as follows. First we 
derive the effective action by the strong couling expan- 
sion. Next we calculate the effective potential (the free 
energy per unit volume at zero temperature) with the 
use of the Hubbard-Stratonovich transformation and the 
mean-held approximation. Finally we obtain the value 
of the chiral condensate as the stationary point of the 
effective potential. 

The main purpose of this study is devided into three 
parts: (I) obtain the value of the chiral condensate, which 
corresponds to the excitation energy of the system, in 
the strong coupling limit. (II) search for the existence 
of the gapless (semimetal) phase between the topologi- 
cal and normal insulator phases. This is due to the fact 
that the transition from topologically nontrivial phase 
to trivial phase (and vice versa) needs the gap closing 
of the bulk energy band. (Ill) search for the phase in 
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which time-reversal and inversion symmetries are sponta- 
neously broken due to electron correlation. Such a phase 
has been confirmed in the lattice quantum chromody- 
namics (QCD) with Wilson fermions 39-41 and was sug- 
gested recently in a mean-field study of Wilson fermions 
with the short-range interaction 42 . 



II. MODEL 

It is known that the effective Hamiltonian of 3D topo- 
logical insulators such as Bi 2 Se 3 is described by the Wil- 
son fermion 5 : 

^o(fc) = ^ . sin kj ■ ctj + m(k)P, (1) 

where m(k) = m + rJ2j (1 — cosfcj), r > 0, j (= 1, 2, 3) 
denotes spacial axis, and ay, /3 are the Dirac gamma 
matrices given by 



er, 

(Ti 



1 

-1 



(2) 



The energy of this system is measured in units of v-p/a 
with vf and a being the Fermi velocity and the lat- 
tice constant, respectively. The Hamiltonian (1) has 
time-reversal (7~) symmetry and inversion (I) symmetry, 
i.e., TUo{k)T- 1 = Ho(-k) and X^o^X" 1 = H (-k) 
are satisfied, where T — 1 ® (—i(J2)fc (/C is the com- 
plex conjugation operator) and X = as <g> 1. In the 
Hamiltonian (1), the spinor is written in the basis of 



-feBf 



'feBJ. 



where is the creation oper- 
ator of an electron, A, B denote two orbitals, and t (4) 
denotes up- (down-) spin 5 . 

In the presence of time-reversal symmetry and inver- 
sion symmetry, the Z 2 invariant of the system is given 
by 3 - 4 



(-ly^IIi-sgnMA,)]}, 



(3) 



i=l 



where Aj are the eight time-reversal invariant momenta. 
It is easily shown that if > too > — 2r or — 4r > too > 
— 6r (mo > 0, — 2r > to > — 4r, or — 6r > to ), the 
system is topologically nontrivial (trivial). 

Let us consider a strongly correlated topological insu- 
lator in Euclidean spacetime, which is described by the 
Wilson fermion with 1/r Coulomb interaction between 
the bulk electrons. We start from the Euclidean action 
of (3+l)D Wilson fermion interacting with electromag- 
netic field on a lattice, which is given by 



+ (to + 4r) ^ i>n4>n, 



(4) 



where P^- = (r ± 7 M )/2. Here n = (no,ni,n2,ns) de- 
notes a site on a spacetime lattice and jj, {ji = 0,1,2,3) 



denotes the unit vector along /x-direction. U n ^ is the 
link variable, which is defined by f7„ ijU = e M 9- 4 M("+/ i / 2 ) j 
where = (A , A) is the four- vector potential, a is 
the lattice constant, and g 2 = e 2 /e with e and e being 
electric charge and the permittivity of the system, re- 
spectively. The term proportional to r is introduced to 
eliminate fermion doublcrs. In this paper, according to 
the Hamiltonian (1), we adopt the Dirac representation 
in the Euclidean spacetime ({7^,7^} = 2S^ U ): 



(5) 



where j = 1,2,3 and aj are the Pauli matrices. 

In the case of 3D topological insulators, the Fermi ve- 
locity vf is about 3 x 10 _3 c where c is the speed of light 
in vacuum. Then the interactions between the bulk elec- 
trons can be regarded as only the instantaneous Coulomb 
interaction (A, = 0) like in the case of graphene 30-38 , so 
the action (4) is rewritten as 

S F = S { F T) + 4 S) + (too + 4r) E V^n, 





'1 " 




—i&j 




'0 1" 


7o = 


-1 


. n = 


i<jj 


, 75 = 


1 



(6) 



where 



4 S) = -V 



(7) 

and U n ,o — e t0n (— n < 9 n < tt). The pure U(l) gauge 
action on a lattice is given by 



S G 



n \x>v 



(Un.jiv + Un,fiv) 



(8) 



where (3 — vp/g 2 - The plaquette contribution U n ^ v is 
defined by 



(9) 



where U n j = 1 in our case. The total action on a lattice 
is written as 



S = Sp + Sg- 



(10) 



The dielectric constant e r of Bi 2 Se3 is rather large 43 (e r = 
e/eo « 100). This means that the Coulomb interaction 
between the bulk electrons in Bi 2 Se3 is considered to be 
weak. In fact, the value of /? is approximated as 



_ v F e r A-Kephc ^ 

P — A ' 9 ~ "J) 



(11) 



and we cannot perform the strong coupling expansion in 
Bi2Se3. However, we think it would be important from 
a theorerical viewpoint to examine the strong electron 
correlation effect in Dirac fermion systems which describe 
topologically nontrivial states. 
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III. EFFECTIVE ACTION 

Let us perform the strong coupling expansion. The 
strong coupling expansion has been often used in 
QCD 44 ~ 48 where the coupling between fermions (quarks) 
and gauge fields (gluons) are strong. We can carry out 
the Uq integral by using the SU(iV c ) group integral for- 
mulae: 



JdUl = l, JdUU ab = 0, J 



dUU ab U} 



cd 



1 



SadSbc- 

(12) 

In the following, we derive the effective action S e f{{ip,ip] 
by carrying out the Uq integral: 



V[^,U ] 



-Sf—Sg 



(13) 



First we consider the strong coupling limit {(3 = 0). In 
this case, the timelike partition function is given by 



7 M 

J SCL 



VU Q e 



3 (t) 



Integration with respect to U is carried out to be 



z SCh = CX P 



(14) 



(15) 



Here we have used the fact that the grassmann variables 
ip's and r/S's satisfy ip 2 = tp 2 = 0. We can rewrite this 
term as 

W VaVd^n = -tr [M n P+M n+b P ] , (16) 

where we have defined (M n ) a p = Vv Q Vv^- The sub- 
scripts a and (3 denote the component of spinors. 

Next we evaluate the term of the order of j3. In order to 
evaluate the plaquette contributions from Sg, we use the 
cumulant expansion 48,49 . Let us define an expectation 
value: 

( A ) = —fry I VU A[U ]e- s *\ (17) 

Then using this definition, the full timelikc partition 
function can be expressed as 



-S G 



Z (T) = J VU e- ■ - ~ SCL 
The contribution from Sq is given by 

AS = -\og(e- Sc ) 



Z[ T X (e~ s °) . 



n=l 



nl 



(18) 



(19) 



where (• • • ) c is a cumulant. The correction to the action 
up to O(fi) is given by 



AS = (S G ) C = (S G ) 



(20) 



The expectation value of U n ^ v is evaluated as follows 48 : 
(U n ^v) - J dU nfi U n ^e- s r\ (21) 

where Sp^ is the plaquette-related part of Sp\ We see 
that the terms with (fi, u) = (i,j) become constant and 
find only (fi, v) = (j, 0) terms to survive: 




(U n ,jo) = -tr 



n,jO 



-tr 



V n,j P V n+ Qj P 
V~ P + V + P~ 



(22) 



and 



where we have defined (V*A , 

V n >J'ap 

Finally, substituting Eqs. (15) and (20) to Eq. (18), 
we obtain the effective action up to 0((3): 

Sett =(m + 4r) ^2 ~ ^2 [$ nP j~^n+j + $„+j P j~lPn 

n n,j 

+ ^tr [M n p+M n+6 P-] 



2 



tr 



V n,j P V n+ oj P 



(23) 



IV. EFFECTIVE POTENTIAL AND CHIRAL 
CONDENSATE 

In this section, we derive the effective potential with 
the use of the extended Hubbard-Stratonovich transfor- 
mation (EHS) 35 ~ 37,48 , and then we obtain the value of the 
chiral condensate as the stationary point of the effective 
potential. We apply the EHS to the trace of arbitrary 
two matrices. Introducing two auxiliary fields R and R', 
we obtain 

eK trAB K 



J V[R,R'\ cxp \ -^Y^ 



(Rafi) 2 + {R'aflY 



(24) 



-(A afj + BTp)Rp a - i(A a p - Blp)R'p a 



where n is a positive constant and the superscript T de- 
notes the transpose of a matrix. Two auxiliary fields 
take the saddle point value, R a p = (^A + B T ) fja and 

R' a p = i(A — B T ) /3a , respectively. Defining Q = R+iR', 
Eq. (24) is rewritten as 



e KtlAB oc 



J P[Q,g*]exp| -kJ2 



\Q a p\ 2 - A a pQi3 a — B^ tj3 Q* !3c 



(25) 
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A. Effective Potential in the Strong Coupling Limit 

We consider to decouple the third term in the effec- 
tive action (23) to fermion bilinear form. To do this, 
we set (k,A,B) = (1, M n P+ , -M n+Q p-) in Eq. (25). 
Assuming that = Q (Q*p a = Q a p) and performing 
the mean- field approximation, Q n _ = Q n = const., we 
obtain 



exp | - Y, tr [MnP + M n+() P -] j 
oc J VQ exp < — Yj 



(Qn) a B -i>n{loQ) tpn 



(26) 

where we have used that [M n ~/ a } afj = ^ n , a ^n^(lo)-,p- In 
this case, Q has the saddle point value 



Q a p = 2<[M„P +] - [M„P -]> 



Pa 



(27) 



We obtain the effective action in the strong coupling 
limit expressed by an complex auxiliary field Q: 



J v[il>,$,Q]t 



(28) 



where 



S cS = N S N T Y \Q a p\ 2 + Y ^kM(k, Q)4> k ^ 

a,/3 k 



with 

M =Y ilj sm % + 
- (7oQ) T • 



m 



•■(4-E, 



cos k. 



(30) 



Here N s = V/a 3 and N T = 1/T with V and T be- 
ing the volume and the temperature of the system, re- 
spectively and we have done Fourier transform from 
n = (n , ni,n2, n 3 ) to k = (fc , fci, k 2 ,k 3 ). 

The effective potential at zero temperature per unit 
spacetime volume is given by 



1 



N S N T 



log Z[Q}. 



(31) 



Integration with respect to ip and tp is carried out by the 
formula / D[if>, $\e~^ M ^ = detM Therefore wc need to 
calculate the determinant of M. Here let us assume that 



(M n ) = cre i075 = cr(cos 91 + i sin 6» 75 ) 



= a 



cos v i sm 9 
i sin 9 cos 9 



Then it follows that 



a cos ( 
cr sin i 



(32) 



(33) 



The terms Oipip} and (V>i 7 5^>) decribe the chiral conden- 
sate and the condensate of pseudoscalar mode, respec- 
tively. The Wilson fermion breaks the chiral symmetry 
by itself (the term proportional to r), hence we use the 
value of (i>ip) to determine the system is whether insulat- 
ing or semimetallic. Namely, if Oipip) 7^ 0, the mass gap is 
dynamically generated and thus the system is insulating. 
From Eqs. (27), (30) and (32), we obtain 



M 



m k - 7 
-Oj sin kj 



A B 
C D 



m k - 



(34) 



where we have defined that m k = m +r ^4 — J2j cos • 
After a straightforward calculation, we have 

detM = detA • dct (£> - CA^B) 

E/in 2 fc J+ (- fc -t) 2 +(t 

' (35) 

Finally we arrive at the effective potential in the strong 
coupling limit: 



J"eff(0<r,</>7r) 



J — 7T 



d 3 k 



x log 



(36) 



The values of (f> a and fa are obtained by the stationary 
conditions dF e s{fa, fa)/dfa = dT^fa, fa)/dfa = 0. 

In the chiral limit (r = too = 0), the effective potential 
is a function of only a, reflecting the chiral symmetry 
of the action. Note that this effective potential corre- 
sponds to that of the staggered fermion (SF) model for 
graphene 35 ' 36 except for an additional factor 4 by setting 
r = and changing from (3+l)D to (2+l)D: 



F e tt{fa, fa) = (<f>v, fa)- 



(37) 



This result is reasonable, because the two cases describes 
the same system where the 2 3 = 8 fermion doublers ap- 
pear. 



B. Chiral Condensate in the Strong Coupling Limit 

At first, we found that the value of fa is zero at the 
stationary point for any set of (r, m ). Hence in the fol- 
lowing, we set fa = —a and 0^ = in Eq. (36) to cal- 
culate the value of the chiral condensate a. The term 
iipj^ip is odd under both time-reversal and inversion. 
Therefore, this means that the phase with spontaneously 
broken time-reversal and inversion symmetries docs not 
arise in the strong coupling (electron correlation) limit. 
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FIG. 1. (Color online) (a) mo-dependence and (b) ./In- 
dependence of the chiral condensate a in the strong coupling 
limit 03 = 0). We define M = m + 4r. 



r = 0.5 









A mean-field study of Wilson fermions with the short- 
range interaction from the weak coupling 42 suggests the 
existence of this phase. Such a phase (where parity and 
flavor symmetry are spotaneously broken) has also con- 
firmed in the lattice QCD with Wilson fermions 39-41 . 

The mo- and M-dependence of the chiral condensate 
<7 is shown in Fig. 1. Here we define M = mo + 4r. 
The value of a is expected to be quantitatively correct, 
based on the fact that the result of a strong coupling 
expansion study in graphene 35 ' 36 is in good agreement 
with that of lattice Monte Carlo studies 30-33 . As men- 
tioned above, at /3 — oo (i.e. non- interacting limit), the 
system with > mo > — 2r (mo > 0) is identified as a 
topological (normal) insulator. The chiral condensate is 
dynamically generated mass gap and corresponds to the 
excitation energy of the system. Therefore, the topologi- 
cal insulator phase changes to the normal insulator phase 
in the strong coupling limit. This is consistent with a 
mean- field analysis from the weak coupling 42 . The value 
of the energy gap in the strong coupling limit is estimated 
as a/2 ■ hv F /a ~ 0.14 [eV] at (r,m ) = (1.0, -0.5), when 
hvp = 3 [eV-A] and a = 3 [A] are assumed. From Fig. 
1(b), we see that a approaches zero as M approaches 
zero and that a gets proportional to 1/M as r becomes 
larger. This r-dependence is similar to that in the lattice 



QCD with Wilson fermions 39 , which is given by 
1 



M 
3M 

16 - M 2 



(\M\ > 2) 
(|M|<2). 



(38) 



However, in the region near \M\ = 2, a is a smooth 
function of M in contrast to Eq. (38). 



C. Effective Potential Up to C»(/3) 

Let us evaluate the G(f3) contribution to the effective 
potential. We write the fourth term in the effective action 
(23) as A5*i + AS^. Then we should choose such that 
(re, A, B) = (/3/2,^„+ P+,-K- +0j P -) in Eq. (24): 



AC" 

e 1 oc exp 



E 



EM( glfej 



R^P+ 



-ik, 



(39) 



Similarly, setting (re, A, B) = ((3/2, V'jP^, -V+ +0 jPf), 
we obtain 



e AS2 cxexp 



-§E 



lj f + & 



E[(*; 

n 



-i e 



-ikj rpijT p+ 



T' j P -) 



4>k 



(40) 

These two equations are very complicated. It is natural 
to suppose the following two cases: 



(i) R — R' and T = T' 

(ii) R = -R' and T = —T 1 . 



(41) 



It is shown that the effective potential of the case (ii) is 
smaller than that of the case (i) [see Appendix]. Thus 
we show the calculation in the case (ii) in the following. 
In this case, combining Eqs. (39) and (40), we have 

ASM,R,T\ =/3]T [{K fi f + (Tip) 2 

j { n 

+ X i>k [-i (R ]T P^ + R'Pv) (42) 



+T 3± P+ -T 3 P 



sin kjipk 
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It have been found in the previous subsection that the 
value of (^751/'), the condensate of pseudoscalar mode, 
is zero for any M in the strong coupling limit. Hence, we 
concentrate on examining the /3-dependence of the chiral 
condensate. For this purpose, we assume that Eq. (42) 
takes time-reversal and inversion symmetric form which 
results from the original Lagrangian. That is, we assume 
that (R° T P^ + R J Po) = 07j and (T^ T P+ - T'Pq) = 
bjj where a and b are constants, since sin kj is odd under 
time-reversal. 

To be concrete, we set as follows: 




er., 



(43) 



where 4>j, ujj are constants. It is easily found that R 2 , 
T 1 and T 3 are symmetric, and that R 1 , R 3 and T 2 
are antisymmetric. However, the final result doesn't de- 
pend on whether the auxiliary field is symmetric or an- 
tisymmetric. Substituting Eq. (43) to Eq. (42), we 
see that the correction from the strong coupling limit 
is obtained by replacing sinfe,- in M(k, cr) (Eq. (34)) by 
(1 — j3r<pj + fiujj) sinfcj. This means that the effect of /3 
is equivalent to the modification the Fermi velocity from 
the strong coupling limit. Then we obtain the effective 
potential expressed by auxiliary fields: 



T cfi (a, (j> 3 , Wj-) = cr 2 + 4/3 YJ$ + u j ) 



n d 3 k 



(44) 



log [detM (k,a,4>j,ujj)] . 



Eliminating <j>j and ujj by the stationary conditions 
dF ee (a,4>j,uij)/d(f)j = dF e s(a,(f>j,LJj)/du)j = 0, we ar- 
rive at the effective potential up to 0{(3): 



F eS (a) =a 2 -2 



(l + r 2 )/3 



3 

E 



d 3 k 
d 3 k 



log 



m k -~ 



sin 2 kj 



E 



sin 2 ki 



iTOfe 



J2i sin2 k i 



(45) 

As well as in the strong couling limit, this effective poten- 
tial gives the same value of the chiral condensate as that 
of the staggered fermion (SF) model for graphene 35 ' 36 
by setting r — and changing from (3+l)D to (2+l)D: 



D. /3-dependence of Chiral Condensate 

The /3-dependence of the chiral condensate a is shown 
in Fig. 2. We see that cr is a monotonically decreas- 
ing function of the coupling strength j3. This behavior 
is consistent with a mean-field analysis from the weak 



(a) 



0.5 
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-0.5r 
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r= 1.0- 



0.4 
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P 

FIG. 2. (Color online) /3-dependence of the chiral condensate 
a with (a) mo = — 0.5r and (b) mo = +0.5r. 



coupling 42 . Our result shows that the mass gap remains 
finite, in contrast to the mean-field analysis in which the 
mass gap becomes infinity in the strong coupling limit. 
We see also that as r becomes smaller, the rate of de- 
crease of a becomes notable. Namely, as the original 
mass of doublers becomes smaller, the energy gap of the 
system becomes smaller, as is understood intuitively. 

It is concluded that the gapped phase (normal insula- 
tor phase) is stable in the strong coupling region. This 
contrasts with the result of the strong coupling expan- 
sion in graphene 35,36 . In graphene, the rate of decrease 
of a from /3 = to j3 = 0.5 is about 60% 36 , whereas that 
of our model is about 10% at r = 0.2. The difference be- 
tween these two cases stems from the spatial dimension 
(two and three). This is because the staggered fermion 
model for graphene and the r = limit of our model 
differ only by the dimension. 



V. SUMMARY 

To summarize, we have studied the strong electron cor- 
relation effect in a 3D topological insulator which effec- 
tive Hamiltonian can be described by the Wilson fermion. 
Based on the U(l) lattice gauge theory, we have per- 
formed the strong coupling expansion. We have found 
that the system always has a finite enery gap in the 
strong coupling region. The values of the chiral conden- 
sate which is regarded as the energy gap of the system in 
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the strong coupling limit are expected to be correct quan- 
titatively. The behavior of the chiral condensate in our 
model is similar to that of the lattice QCD with Wilson 
fermions. The phase where time-reversal and inversion 
symmetries are spontaneously broken was not found in 
the strong coupling limit, in contrast to the case of the 
lattice QCD. 

The critical value of coupling strength, at which the 
system changes from insulator to semimetal, could not 
be obtained, although it was shown that the energy gap 
becomes smaller as electron correlation gets weaker. For 
the consistency with the result in a weak coupling anal- 
ysis, it is concluded that the topological insulator phase 
changes to the normal insulator phase in the strong cou- 
pling limit. In this study, the bulk property of a 3D 
topological insulator was examined. It will be interest- 
ing to examine the strong correlation effect in the surface 
Dirac fermions. 
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Appendix A: Determination of Effective Potential 

In this appendix, we show numerically that the effec- 
tive potential of the case (ii) is lower than that of the 
case (i). In the case (i), R = R' and T = T', combining 
Eqs. (39) and (40), we have 

AS[^,R,T] =13 IE [«0 3 + ( T ^) 2 " 

j { n 

k 

+i (T^ T P+ + T^Po) ] cosfcjVfc j ■ 

(Al) 



Since cos kj is even under time-reversal, we may choose 
such that 



(A2) 




where (j)j, ujj are constants. In this case, we see that the 
correction from the strong coupling limit is equivalent to 
the modification of the masses of doublers. That is, the 
effective action is obtained by replacing cos kj in j\4(k, a) 
(Eq. (34)) by (1 + /30j + rfiujj) cos kj. Then we obtain 
the effective potential expressed by auxiliary fields: 



-^off (c, 1 



3 ' > 



4/3E( 



w d 3 k 



(A3) 



log [det M (fc, cr, (j)j , ujj )] 



Eliminating <j)j and ujj by the stationary conditions 
dF eS (a,<f>j,ojj)/d<j>j = dFcsia^j^uj^/dujj = 0, we ar- 
rive at the effective potential up to 0(f3): 



F eS (a) = <j 2 -2 



(l-r 2 )/3E 



K d 3 k 

* d 3 k 



log 



( m fe _ f ) + E sin2 ki 



{mk — f ) cos kj 



+ Hi sin2 k i 



(A4) 

Comparing Eq. (A4) and Eq. (45), we see that the 
ground state is realized in the case (ii) [see Fig. 3]. 
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fen 
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-2.72 
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FIG. 3. (Color online) Comparison of F e g(a) in the case (i) 
and (ii). r = 0.50, m = -0.50, /3 = 0.30. 
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